Bona Fide Thermodynamic Temperature in Nonequilibrium Kinetic Ising Models 
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We show that a nominal temperature can be consistently and uniquely defined everywhere in the 
phase diagram of large classes of nonequilibrium kinetic Ising spin models. In addition, we confirm 
the recent proposal that, at critical points, the large-time "fluctuation-dissipation ratio" Xtx> is a 
universal amplitude ratio and find in particular Xoa « 0.33(2) and X<x> = | for the magnetization 
in, respectively, the two-dimensional Ising and voter universality classes. 
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The fluctuation-dissipation theorem (FDT), which re- 
lates the correlation and response (or susceptibility) func- 
tions during the return to equilibrium following a small 
perturbation, may be used to provide an absolute defini- 
tion of the temperature of a physical system. However, 
there are many instances, such as glassy materials and 
coarsening systems, where the standard form of the FDT 
breaks down, because the equilibration timescales are as- 
tronomical or even infinite in the thermodynamic limit. 
Relaxation properties then depend upon both t, the time 
at which perturbations are applied, and t' , at which mea- 
surements are taken, giving rise to "aging phenomena" . 
In this context, guided by the dynamics of some mean- 
field spin glass models, Cugliandolo and Kurchan pj have 
proposed a generalized form of the FDT: 



R{t,t') = 



X(t,t') dC{t,t') 



T 



dt 



(1) 



where the two-time response and auto-correlation func- 
tions R(t,t') and C(t,t') of some physical observable, as 
well as the so-called FDT ratio X(t, £'), are not functions 
of the t' — i, and, moreover, X depends functionally on 
C . The limit value X^ often taken by X at large times 
is usually interpreted as the emergence of an "effective 
temperature" T oS = T/X^ @. 

The picture suggested by the above generalized FDT 
scenario has received a large number of confirmations Q , 
both from experimental 4] and numerical |5|,|6|)l3 studies 
of structural and spin glasses, granular matter, and gen- 
tly sheared fluids, to quote a few examples. In these situ- 
ations, the physical systems under consideration, though 
they may never achieve equilibration, evolve subject to 
a temperature T. However, this "nominal" temperature 
is generally not even defined for intrinsically nonequilib- 
rium systems, e.g. stochastic models for which detailed 
balance is violated, although such systems are known 
to behave much like their equilibrium counterparts, in- 
cluding critical behavior and non-stationary properties of 
their two-time autocorrelation and response functions. 

In this Letter, we show that a nominal temperature 
can be consistently and uniquely defined everywhere in 
the phase diagram of large classes of nonequilibrium ki- 



netic Ising models. We explain how both a "dynamical" 
definition of temperature (using Eq. Q) and a "geomet- 
rical" approach (in terms of the density of states) lead 
to an ambiguity that can be lifted using a maximum- 
entropy argument. This results in a unique, bona fide 
thermodynamic temperature which takes the same value 
for different observables and coincides with the usual one 
when detailed balance is satisfied. In addition, consid- 
ering the critical points of our models, we confirm the 
proposal 0, E| that Xqo is, in this context, a universal 
amplitude ratio of dynamic scaling functions and find in 
particular I M [M] w 0.33(2) and X^M] = \ for the 
magnetization M in, respectively, the two-dimensional 
Ising and voter classes. 

For simplicity, we illustrate our approach by con- 
sidering the family of two-dimensional nonequilibrium 
kinetic spin models introduced some time ago in [llj 
which are defined by the following evolution rules. Dur- 
ing an elementary timestep, an Ising-like spin (oy = 
±1) on a square lattice is randomly picked up, and 
flipped with a probability W[E r = —^a r H r ] where 
E r € {—2, —1, 0, 1, 2} is the local (pseudo-) energy, H r . = 
^] j Ur+e^ being the local (Weiss-like) field calculated 
over the four nearest neighbors of oy. The local Z2- 
symmetry of the dynamical rules is enforced by demand- 
ing that W{— E) = 1 — W(E) (the system being homo- 
geneous, we drop spatial indexes whenever possible) and 
thus W(0) = \, leaving a two-parameter family defined 
byH/(l) = \{l + x) andl^(2) = ±(l + y), with —to favor 
ferromagnetic ordering — < x , y < 1. An alternative 
convenient parameterization 12] is to match the flip rate 
W with the corresponding expression for Glauber dy- 
namics (W = |(l + tanh2/3_B) at temperature T = 1/(3). 
This is only possible if one introduces two parameters 
Pi = (3e=i and @2 = Pe=2, such that x — tanh2/3i and 
y = tanh4/?2, which measure respectively the strength 
of interfacial and bulk noise 0|. The usual Glauber 
dynamics then corresponds to (3\ — fi% — l/T, or, equiv- 
alently, y — 2x/(l + x 2 ). Other well-studied models in 
the (x, y)-plane comprise the majority model {y = x), 
the "noisy voter" or linear model (y = 2x), and the "ex- 
treme" model (x = 1). Except on the Glauber line, mod- 
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els in the (x, y)-plane do not obey detailed balance or pos- 
sess an underlying short-range Hamiltonian. Neverthe- 
less, in agreement with a long-lastin g co njecture backed 
up by field-theoretic arguments |l5l llq . there exists a 
critical line y c (x) separating a disordered (paramagnetic) 
from an ordered (ferromagnetic) phase and along which 
Ising static and dynamic critical exponents are numeri- 
cally found (Fig. Hi) HHHE!. This line terminates 
at the voter model critical point (|, 1), across which the 
transition occurs in the absence of bulk noise (y = 1) 

The gist of our argument to define an effective tem- 
perature is the observation that, even if the FDT may be 
broken at large times, there always exists, for small time- 
differences (if — t < 0(1) and t 3> 1), an equilibrium- 
like regime where the standard form Q is valid with 
X = 1. Thus, in this two-time sector, and in particular 
for t' — t <C 1, we view as a means of obtaining a 
dynamical definition of an effective temperature |l7| : 



(a) 



T, 



1_ = lm R(tt>) 

t->oo,i'-i->0 dtC(t,t') 



dyn 



(2) 



To really make sense of {2J , we must of course specify the 
response function R(t,t'), i.e. we have to implement, for 
a generic (x, y)-model, the effect of a symmetry-breaking 
external field h r linearly conjugated to the local spin vari- 
able ay. Probably the simplest choice is to make the 
substitution a r H r — > a r (H r + h r ) in the original expres- 
sion of W r , this simple recipe being also consistent with 
the natural redefinition E r — » E r ^ ' = E r — a r h r of the 
energy in this case. This yields: 

W r (h) = \ (1 - o- r tanh [P\ Er \(H r + h r )}) . (3) 

The problem we face now is that an extra parameter /3 
appears, which governs the fate, under the influence of 
the applied field h r , of configurations where H r = 0. One 
might think that the arbitrariness in the choice of (3q for a 
generic (x, y)-model could be resolved by demanding that 
(0 would give the same effective temperature for all ob- 
servables. This is not the case: in fact, general properties 
of the spin-flip dynamics we consider automatically guar- 
antee that, for any given value of fto, our definition (J2J of 
Pdyn = 1/Idyn does not depend on the observable chosen 
to evaluate the corresponding correlation and response 
functions since it can be re-expressed as 



dyn 



(a r (d h Wy i) )h=o) _ (2(3 lErl W r (l-W r )) 



(W r ) 



(W r ) 



(4) 



where the average is over all lattice sites . 

The way out of this quandary is deceptively simple: as 
for the Glauber model, where @o is of course equal to the 
inverse of the true temperature, one has to tune /?o such 
that any FDT measurement with a probing field yields 
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FIG. 1: (a): phase diagram in the (x, y)-plane: the Ising- 
critical line separating the disordered region (below) from the 
ordered one (above) terminates at the voter critical point (big 
filled circle), (b): /3d yn (crosses) and /3 goo (diamonds) along 
the linear model line, (c) Inverse effective temperature (com- 
puted from /3dyn, exact for Glauber and voter) along the crit- 
ical line. 



precisely /?o as the effective temperature. That is to say, 
for any (x, y)-model, /3q is the solution of 



/?dyn(/3o) — 00- 



(5) 



Because /?□ only appears, in Eqs.©-®, in the initial 
condition R(t,t + ) of the response function ^^|, the de- 
pendence of /?dyn on (3 is simply linear, and the solution 
of (J3J) is therefore unique. 

The reason why Eq.([5j is the correct prescription is a 
profound, geometric one, and is incidentally the same as 
for the Glauber model: this choice of (3q actually maxi- 
mizes the entropy of the spin configurations. To under- 
stand this point, recall that at equilibrium the canon- 
ical Botzmann-Gibbs distribution can also be obtained 
from the microcanonical ensemble under an information- 
theoretic, maximum-entropy condition |2(ij |. where the 
Lagrange multiplier (3 gco (identified with the inverse tem- 
perature in units where kg — 1) imposing the mean value 
of the energy reweights with an exponential prefactor 
cx e /3gcoE the bare probabilities of the spin configura- 
tions. Then, for any (x, y)-model on a finite system a 
steady-state is always attained, and it is perfectly legit- 
imate to define a microcanonical entropy S m ic (M, E) = 
lnO(M, E) where the density of states fi(M, E) simply 
counts the number of configurations having a prescribed 
magnetization M and energy E as reached by the spin- 
flip dynamics. But in this "microcanonical" ensemble, 
the only spin-flip allowed are those which do not change 
the energy, that is the ones corresponding to W(Q) = \, 
a value common to all (x, y)-models. In other words, the 



influence of distinct (x, y)-values can only be felt if one 
lets the energy fluctuates, in which case — and for ex- 
actly the same reasons as above — one should expect 
to actually observe the configurations with an effective 
reweighting oc e^ scoE , leading to a "canonical" entropy 
Scan- Now, if one works with a non-zero external field 
h in the microcanonical ensemble (in particular around 
the most probable values of M and E), to lowest-order 
in h, once again, only the spin-flips with H r = will be 
allowed, with a probability © which can be rewritten 
as W^(0) = e it o h^-0 o h ■ This expression shows that 
a consistent reweighting can be achieved if and only if 
figco — Po, an equality that (JjjJ also recovers if one re- 
stricts the average to zero local-field configurations. 

To check the validity of the above scenario, we have 
computed the density of states (reweighted) using an ef- 
ficient algorithm |2JJ which allows to determine, with 
high accuracy, the ratios q^'^) eP°^ E ~ E - > for neighbor- 
ing (i.e. connected by a single spin-flip) magnetization 
and energy levels. This method is applicable anywhere in 
the (x, y)-plane, except of course for y = 1, 1/2 < x < 1 
when no fluctuations are present. It produces, in partic- 
ular, fast and accurate estimates by restricting the calcu- 
lations to a narrow interval centered around the station- 
ary (and most probable) values M sta and i? s ta, with or 
without an external field. The above ratios also provide 
a direct access to an estimate of as ^(M^f3 h) wmc j 1 _ 
if the above effective thermodynamic picture is valid — 
should be equal to h/3 gco in the disordered phase (where 
for any system size Af sta = 0). This is indeed verified for 
sufficiently small h and, most importantly, it is only so 
when one has tuned (3q according to (Fig.[2t). 

In the ordered phase or along the critical line, finite- 
size fluctuations of the energy bring in another term 
(stemming from the differentiation of e^ sc ° Eih) in the ef- 
fective occupation probabilities) which cannot be evalu- 
ated simply with the present method. But one can nev- 
ertheless have an indirect access to the dependence of 
the absolute entropy on (3q by working, with and with- 
out magnetic field, along the magnetization direction in 
an interval comprising both M^l and , because this 
last quantity depends on (3q- In particular, if the se- 
lected value of (3q is the extremalizing one, the difference 
SS(M) = ^ME, M _ as c „(M, g ,o) should be as flat 

as possible. This is confirmed in Fig. {2fc> for the critical 
point of the majority model. 

To sum up, our nominal temperature is unique and 
can be determined in two equivalent ways for any (x, y) 
model (see Fig. ^ f° r an example). We note in pass- 
ing that two models sharing the same stationary energy 
(e.g. Glauber and linear for which this quantity is known 
exactly H3]) are usually not at the same effective tem- 
perature. Of particular interest is the non-monotonous 
variation of the temperature T c along the critical line 
(Fig. 2J;) . Once T c determined, one can study the long- 
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FIG. 2: (a): for the linear model at i = 0.1 (disordered 
phase), dS ca .n(M, h(3o)/dM evaluated at M sta = for various 
h values for /3 = 0.02, 0.0501, and 0.08 (from bottom up). It 
is only equal to h/3o (solid lines) for the middle value, equal to 
/3dyn- (b) 5S vs M for the critical majority model (x c ~ 0.850) 
with p = 0.365, 0.415, and 0.465 from bottom up (h = 0.1, 
M sta — 0.666, dotted lines are linear fits). Stationarity is 
achieved for the middle value, equal to Pd yn - 



time behavior of X. Using a now-standard method [23j, 
the parametric representation of Fig|3^ yields the univer- 
sal limit Xoo = 0.33(2). This value while in disagreement 
with [tj , has also been obtained independently in 1241 l25l , 
and is consistent with the two- loop RG results of 26]. 

Further (analytical) insight can be gained for the lin- 
ear and critical voter model: along the line y — 2x, all 
correlation functions do not couple to higher-order ones, 
and verify diffusion equations. Using standard Laplace- 
Fourier methods much as for the Glauber ld-Ising chain, 
one can not only calculate exactly the stationary energy, 
but also obtain many exact results for the correlation and 
response functions. Skipping all technical details [22| . 
the two-time two-point correlation function C r (t,t') = 
(ao(t)a r (t')), is compactly expressed in terms of a dou- 
ble Laplace transform C r (s,s') = - } /n \ ^t^/ 2 )-^-^ ) 

^ ' V ' / sqo(s/2) s'~s/2 

(s <->• t, s' <-> t'), where q r (t) = e- {1 - 2x ^p r {2xt), p r (t) be- 
ing the probability that a simple random walk goes from 
to x during a length of time t (po{t) — e~ l 7q(£/2) ~ 
l/(7rf), where Iq is a modified Bessel function). Special- 
izing to the autocorrelation function C(t,t') = Co(t,t'), 
this gives for x < 1/2 when t is large and t' — t arbitrary 

C(t,t') ~ 1 — f Q 1 duq (u) j qo(0) , a function of t' — t de- 
caying exponentially on a timescale t x = 1/(1 — 2x) . At 
the voter critical point, both <?o(0) and t x diverge, and the 
previous expression crosses over to a non-trivial two-time 
scaling form, with C(t,t') ~ 1 — (hit) 1 f Q * dupo(u) in 
the first temporal regime t 3> l,t' — t ~ 0(1) , which 
matches at large time differences the result C(t,t') = 
(kit) - In (frrf) valid in the second regime t,t'-i>l. 
As for the autoresponse function R(t,t') it simply reads 
R{t,t') = q (t' - t)R(t,t+). For x < 1/2, R(t,t + ) 
tends to a finite and non-zero constant, and an effec- 
tive FDT is obeyed at all times, with X(t,t') = 1. For 
the voter model, R{t,t+) ~ (tt / lnt)[f3 - ^(/3 - Pi)} 
in the long-time limit, where the numerical constant C4 
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FIG. 3: (a) Parametric plot of x/(C/3) vs. C for the Glauber, 
majority, and extreme models at their critical point, where 
C = C(t, f) is the spin autocorrelation function, and x(t, t') = 
f duR(u,t') is the integrated response following the applica- 
tion, on random initial conditions, of a bimodal random field 
of magnitude h m 0.05 switched off at time t = 50 (single 
runs on lattices of 8192 2 spins with periodic boundaries, data 
shown for t' € [t, 3000] For small C, all curves level off around 
X X [M] = 0.33(2). (b) Short-time part of the x - C plot for 
the local magnetization (bottom) and the local energy (top) 
in the voter model. Both curves have slope /3 V otcr = ^p- 

is such that (n =1 er M (t)) - 1 - c 4 7r/lnt. Hence |@} 

reads /3dyn(/3o) = Po ~ §(A) — Pi)c4, and the prescrip- 
tion (jSJ is realized when (3q = thus bypassing the 
evaluation of C4. Therefore the strength of interfacial 
noise driving the voter order-disorder transition can be 
unambiguously associated with a temperature exactly 
given by l/T c = ^ = 0.27465..., and the FDT ra- 
tio X(t,t>) * po[ /4* + -%„ t) ^X ao = \ 21\. While T c 
is easily measured for both M and E and found equal 
to the exact value above (Fig. Et>), the long-time regime 
is difficult to reach numerically. For the magnetization, 
we find X^M] = 0.50(5) (not shown), which gives an 
indication of the accuracy of the Ising value above. For 
the energy, however, the signal is too weak to estimate 
Xoo[E] so that we cannot test the conjecture that 
is the same for all observables, which seems unlikely 
to us given that critical static amplitudes ratios usually 
depend on the observable considered. 

Summarizing, we have shown how to measure the 
"strength of noise" via a bona fide thermodynamic tem- 
perature in kinetic spin models and confirmed the univer- 
sality of the FDT ratio at critical points. Ongoing 
work aims at extending our method to even more general 
nonequilibrium systems such as chaotic map lattices. 
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